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Abstract: Porcine deltacoronavirus (5-CoV) is the object of extensive research in several countries 
including the United States. In contrast, the epidemiology of 5-CoVs in wild birds in the US is largely 
unknown. Our aim was to comparatively assess the prevalence of 5- and y-CoVs in wild migratory 
terrestrial and aquatic birds in Arkansas, Illinois, Indiana, Maryland, Mississippi, Missouri, Ohio, 
Tennessee and Wisconsin. A total of 1236 cloacal/fecal swabs collected during the period 2015-2018 
were tested for y- and 5-CoVs using genus-specific reverse transcription-PCR assays. A total of 61 
(4.99%) samples were y-CoV positive, with up to 29 positive samples per state. In contrast, only 14 
samples were positive for 5-CoV (1.14%) with only 1-4 originating from the same state. Thus, unlike 
previous reports from Asia, y-CoVs are more prevalent than 5-CoVs in the US, suggesting that 
5-CoVs may spread in birds with lower efficiency. This may indicate 5-CoV emerging status and 
incomplete adaptation to new host species limiting its spread. Phylogenetic analysis of the partial 
N gene revealed that the newly identified 5-CoV strains were most closely related to the HKU20 
(wigeon) strain. Further studies are necessary to investigate the role of aquatic bird 5-CoVs in the 
epidemiology of 5-CoVs in swine and terrestrial birds. 


Keywords: coronaviruses, 5-coronavirus; wild birds; epidemiology; United States 


1. Introduction 


Coronaviruses (CoVs) are positive-sense RNA viruses that are widespread in humans alongside 
various mammalian and avian species. They cause enteric, respiratory, or systemic diseases of variable 
severity [1,2]. CoVs belong to the Coronaviridae family that contains four genera: Alphacoronavirus 
(x-CoV), Betacoronavirus (B-CoV), Gammacoronavirus (y-CoV) and Deltacoronavirus (5-CoV) [3]. The virus 
evolves through an accumulation of point mutations and both homologous and non-homologous 
recombination. It is hypothesized that this ability to recombine allows the virus to evolve and 
create new forms which can target different species [4]. ~-CoV and £-CoV infect multiple species 
of mammals, but y-CoV is only known to infect birds. In contrast, 5-CoVs are found in both birds 
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and mammals including pigs [4]. Porcine deltacoronavirus (PDCoV) emerged in the US in 2014 as a 
new enteropathogenic CoV causing diarrhea, vomiting, and mortality in neonatal piglets, resulting in 
economic losses to the swine industry [4,5]. Its origin is unknown, but it is hypothesized that it may 
have spilled over from an avian host. 

5-CoV and y-CoV are co-circulating in wild avian species, but their epidemiology (prevalence and 
diversity) differ [6-13]. In a screening study of 918 wild Australian birds, 141 of them tested positive 
for CoVs. After sequencing, 5-CoV was detected in pacific black ducks, curlew sandpiper, red-necked 
stints, ruddy turnstones, and pied herons [6]. However, y-CoV (80% of sequenced samples) was 
identified more frequently than 5-CoV (20% of sequenced samples) [6]. y-CoV was also detected in 
quails and pheasants in Italy [7]. RNA-dependent RNA polymerase (RdRp) analyses found that quail 
were also infected with 5-CoV [7]. Wild birds harbor genetically diverse 5-CoV with some of them 
potentially transmittable to pigs [8]. A previous report from Brazil showed that 5-CoVs were detected 
in two species of wild terrestrial birds, purple-breasted-parrot and plain parakeet. These avian CoVs 
were monophyletically related to CoVs from Sparrow (SPaCoV HKU19) and swine (PorCoV HKU15), 
which were also found to be 5-CoV. This information supports the hypothesis that 5-CoV can cross 
inter-species barriers, with the potential to transmit from birds to swine [8]. 

Previous studies of CoVs in Hong Kong reported three novel CoVs in bulbuls (BuCoV HKU11), 
thrushes (ThCoV HKU12) and munias (MuCoV HKU13), which were hypothesized to belong to the 
novel genus 5-CoV [9]. These three novel 5-CoVs clustered with a 5-CoV detected in Asian leopard 
cats [10]; however, further studies are necessary to understand interspecies transmission from birds [9]. 
Ina surveillance study in Hong Kong and Cambodia, 5-CoVs were found in different wild aquatic birds: 
gray herons, pond herons, great cormorants, black-faced spoonbills, and several duck species [11]. 
5-CoVs do not cause severe illness in birds, leading to their endemic nature in the avian population [11]. 
In a screening study for 5-CoVs in China, none were found in Asian leopard cats, bats, domestic 
cats, cattle, chickens, dogs, humans, monkeys, and rodents [12]. The investigators discovered seven 
novel 5-CoVs in pigs and wild birds (white-eye, sparrow, magpie robin, night heron, wigeon, and 
common moorhen) [12]. Genome sequencing and comparative genome analysis showed that avian 
and swine 5-CoVs had similar genome characteristics and structures [12]. Similarly, a recent study in 
the US described a novel sparrow 5-CoV that clustered together with PDCoV [13] and other terrestrial, 
but aquatic bird 5-CoVs. 

In Brazil, screening revealed the presence of CoVs in vinaceous-breasted amazon and plain 
parakeet, that were closely related to 5-CoVs from birds (SpaCoV HKU17) and swine (PorCoV). This 
emphasizes the increased risk for direct interspecies transmission and that in contrast to waterfowl, 
terrestrial birds may act as intermediate hosts and are commonly found in the sites of bird-swine 
co-mingling [8]. In Finland, two different avian 5-CoV strains from lesser black-backed gull and 
black-headed gull were detected which shared 83% and 85% nucleotide identity with avian and 
mammalian 5-CoVs, respectively [14]. 

In summary, recently identified 5-CoVs are emerging globally that possess high genetic and 
antigenic plasticity. Porcine 5-CoV is being studied extensively in the US and elsewhere [4,5]. In contrast, 
the diversity and epidemiology of 5-CoVs in wild birds in the US are largely unknown. The association 
of porcine 5-CoV with enteric disease in pigs and its continuous spread underscore an urgent need to 
further investigate the mechanisms of its persistence and spread. Extensive epidemiological studies 
are needed to evaluate 5-CoV prevalence in avian species, swine farms, pig—bird comingling sites and 
high pig-traffic sites such as abattoirs, assembly yards and truck washes. In this study, we investigated 
the prevalence of ô- and y-CoVs in wild migratory terrestrial and aquatic birds in Arkansas, Illinois, 
Indiana, Maryland, Mississippi, Missouri, Ohio, Tennessee, and Wisconsin. 
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2. Materials and Methods 


2.1. Samples and RNA Extraction 


A subset of 1236 avian cloacal/fecal swabs were selected from a collection of 16,672 that were 
collected for influenza A virus surveillance in different species of wild terrestrial and aquatic birds 
during the period 2015-2018. A total of 1236 avian cloacal swabs were collected for avian influenza virus 
(AIV) surveillance in different species of wild terrestrial and aquatic birds along the Mississippi Flyway, 
US during the period 2015-2018. Of 1236 avian cloacal swabs, 736 were from aquatic bird samples 
collected in Ohio, Illinois, Missouri, Maryland, Wisconsin and Tennessee in the period 2017-2018. 
A total of 500 were from terrestrial and aquatic bird samples collected in Ohio, Mississippi, Indiana 
and Arkansas in the period 2015-2016. 404 were from known species: 234 aquatic and 170 terrestrials; 
the remainder were from environmental fecal samples. RNA from the avian cloacal/fecal swabs was 
extracted using a modified commercial protocol (Ambion® MagMAX™, Applied Biosystems, Foster 
City, CA, USA) with 50 mg/mL of bovine serum albumin (BSA) and 17% sodium sulfite. 


2.2. y- and 6-CoV RT-PCR Assays 


A total of 1236 avian cloacal/fecal swab RNA samples were tested for y- and ô-CoV 
using pancoronavirus (IN2deg/IN4deg) and deltocoronavirus-specific reverse transcription-PCR 
(RT-PCR) assays, respectively [15,16]. 5-CoV primers were designed using ô-CoV 
nucleocapsid (N) universal primers (UDCoVF: 5’-RYWGAYKSNTCNTGGTTYCA-3’ and UDCoVR: 
5’-HGTGCCWGTRTARTARAAGG-3’) targeting 194bp [16]. Subsequently, after evaluating the results 
of Next-Generation Sequencing (NGS), we designed HKU20-specific primers targeting a 384 bp 
fragment of the N gene (HKU20-N-F24314 5’-TCCGCGCCTCATGGCTCTC-3’ and HKU20-N-R24698 
5’-TCATGAGAAGGATTCTAG-3’). All the RT-PCR assays were performed under the same conditions 
using QIAgen one-step RT-PCR kit (QIAGEN Inc., Valencia, CA, USA) in a GeneAmp PCR system 
9600 thermal cycler (Applied Biosystems, Foster City, CA, USA). The reaction-mixture (total 25 uL) 
included 5x QIAGEN OneStep buffer (5 uL), dNTP (1 uL), upstream and downstream primers (100 
umol/L, each 1 uL), RNAsin (40 U/uL, 0.25 uL), enzyme mix (1 uL), RNase-free water (10.75 uL) and 
RNA template (5 uL). The thermocycler protocol consisted of an initial reverse transcription step at 
50 °C for 30 min, followed by PCR activation at 95 °C for 15 min, 40 cycles of amplification (95 °C 
for 20 s, 50 °C for 10 s, 60 °C for 30 s), and a final extension step at 72 °C for 10 min. PCR products 
were analyzed on 2% agarose gel. The PCR product was purified with an agarose gel QIAquick Gel 
extraction kit (QIAGEN Inc., Valencia, CA, USA). 


2.3. Sequencing and Sequence Analysis 


Amplicons derived using the HKU20-specific primers were sequenced at the Genomics Shared 
Resource of The Ohio State University (Columbus, OH, USA) by the Sanger method [17]. For 
NGS, previously extracted RNA underwent cDNA synthesis according to a random primer protocol 
performed using RevertAid H Minus First Strand cDNA Synthesis Kit (Thermo scientific, Waltham, 
MA, USA). PCR was conducted using True-Start DNA polymerase with 10 mM dNTPs mix and 10 
pmol specific primers per reaction (Thermo Scientific), according to the manufacturer’s protocols. 
TruSeq Stranded Total RNA Library Prep Kit was used with 1 ug total RNA for the construction of 
libraries according to the manufacturer’s protocol. For rRNA-depleted library, rRNA was removed 
from 2.5 ug total RNA using Ribo-Zero rRNA Removal Kit (mixture 1:1 Human/Mouse/Rat probe and 
Bacteria probe), according to the manufacturer’s protocol (with probe concentration for epidemiology 
kit protocol). All cDNA libraries were sequenced using an Illumina HiSeq2000 (Illumina, San Diego, 
CA, USA), producing 101 x 7 x 101 bp paired end reads with multiplexing. Reads were trimmed 
using default parameters with CLC Genomics Workbench 8.5.1 (Qiagen Bioinformatics, Redwood City, 
CA, USA). Trimmed reads were de novo assembled using a word size of 64, bubble size of 100, and 
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minimum contig length of 100. The contigs were subject to the BLASTN search. 5-CoV sequences were 
deposited into GenBank with the accession numbers MN379902, MN379903 and MN379904. 

Sequences were assembled using BioEdit Sequence Alignment Editor and aligned using ClustalW. 
The N gene sequences of avian and porcine 5-CoVs obtained through BLAST search and from GenBank 
were included in the phylogenetic analysis. A maximum likelihood phylogenetic tree was constructed 
using MEGA7 software [18]. 


3. Results 


3.1. Detection of y-CoV in Wild Birds 


A total of 61 out of 1236 birds screened were positive for y-CoV with an apparent prevalence of 
4.99% (95% confidence interval: 3.9%-6.35%). y-CoV was identified in the states of Missouri (n = 29), 
Wisconsin (n = 10), Illinois (n = 8), Tennessee (n = 7), Ohio (n = 5), and Maryland (n = 2) (Figure 1A). 
The states where y-CoVs were identified are shown in Figure 1A. y-CoVs were detected in 6 different 
bird species: blue winged teal (Spatula discors) (n = 27), mallard (Anas platyrhynchos) (n = 16), American 
green-winged teal (Anas crecca) (n = 15), northern pintail (Anas acuta) (n = 1), ring-necked duck (Aythya 
collaris) (n = 1), and American wigeon (Mareca americana) (n = 1), all of which are aquatic bird species, 
meaning a prevalence of 6.3% in aquatic birds and 0% in terrestrial birds. Of interest, y-CoVs were 
detected most frequently in samples from American green-winged teal, blue-winged teal and mallard 
that were among most frequently sampled waterfowl species (Table 1). 


Latitude 
Latitude 


-50 E 
Longitude Longitude 


Figure 1. Gammacoronavirus (A) and Deltacoronavirus (B) detection in different states of the Mississippi 
Flyway (blue contour) and Atlantic Flyway (red contour). 


Table 1. Frequency of y-CoV and ô-CoV in individual bird species. 


Bird Species y-CoV+ 6-CoV+ Total 
American green-winged teal 27/24.5% 0/0% 110 
American wigeon 1/5.6% 0/0% 18 
Blue-winged teal 27/21.4% 6/4.8% 126 
Mallard 16/7% 2/0.9% 227 
Northern pintail 1/8.3% 0/0% 12 
Northern shoveler 0/0% 1/16.7% 6 
Red-tailed hawk 0/0% 1/100% 1 
Ring-necked duck 1/4.3% 0/0% 23 


Snow goose 0/0% 4/5.8% 69 
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3.2. Detection of 6-CoV in Wild Birds 


There were 14 of 1236 birds that screened positive for 5-CoV, corresponding with an apparent 
prevalence of 1.13% (95% confidence interval: 0.68%-—1.91%). Thus, unlike previous reports from 
different countries [6,11,14], our study showed that in the US, y-CoVs are more prevalent than 5-CoV. 
5-CoVs were detected in the states of Illinois (n = 4), Arkansas (n = 4), Ohio (n = 2), Missouri (n = 2), 
Maryland (n = 1) and Mississippi (n = 1) (Figure 1B). The positive sampling locations are shown in 
Figure 2. The bird species identified were the blue winged teal (Spatula discors) (n = 6), snow goose 
(Anser caerulescens) (n = 4), mallard (Anas platyrhynchos) (n = 2), red-tailed hawk (Anser caerulescens) 
(n = 1) and northern shoveler (Spatula clypeata) (n = 1). Of these 14 positive samples, only up to 4 
originated from the same state (Figure 1B), suggesting that 5-CoVs spread with low efficiency in the 
avian species tested. Of note, the red-tailed hawk was sampled upon intake to a wildlife rehabilitation 
center. Interestingly, the prevalence of 5-CoV in aquatic birds was 1.34% compared to only 0.6% 
in terrestrial birds, suggesting that aquatic birds represent an important natural reservoir for CoVs. 
Blue-winged teal and snow goose were the two species associated with the most frequent recovery of 
ô-CoVs (Table 1). 
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Figure 2. Phylogenetic trees based on partial sequences (generated using NGS) of the N (A), S (B) and 
M (C) genes of Blue Winged Teal coronavirus/USA/Ilinois2562/2017. 


3.3. Sequence and Phylogenetic Analyses of 6-CoVs 


In contrast to y-CoVs that are widespread and have been circulating in the US for decades, the 
prevalence and genetic characteristics of ’-CoVs in the US are not known. To address this gap, we 
conducted sequence and phylogenetic analysis of the 5-CoVs identified in this study. All amplicons 
generated with degenerate UDCoV-specific primers were confirmed to contain -CoV N-gene sequences 
sharing ~69%-80% nucleotide identity with various avian 5-CoV species. However, due to a high 
number of degenerate nucleotides incorporated in these primers to allow for broad reactivity and 
detection of highly diverse porcine and avian 5-CoVs, the amplicons were not of satisfactory quality, 
with up to 10% of ambiguous nucleotides (Ns). Thus, we selected 4 5-CoV positive samples (Table 2) 
that had sufficient amounts of RNA and subjected them to NGS sequencing (Table 2). For one of the 
samples, Blue Winged Teal coronavirus/USA/Ilinois2562/2017, NGS recovered approximately 42% 
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of eukaryotic (host genome; normal; incidental), 20% bacteria, 5% virus, and 33% other. Bacterial 
reads identified Fusobacterium nucleatum (18%), while further analyses of viral and “other” reads 
identified various bacteriophage (80%). Eighty-five reads were identified as 5-CoV including fragments 
of ORFla/b, S, M, N and NS7a genes. Blast search and phylogenetic analysis identified that these 
genomic fragments shared the highest nucleotide identity (87%-93%) with Wigeon CoV HKU20, 
without evidence of recent recombination events (Figure 2A—C). For the other three samples, there were 
variable amounts of eukaryotic (2%-74%), bacterial (1%-20%), viral (5%-92%) and other/bacteriophage 
(5%-33%) reads; but no 5-CoV sequence was recovered. This suggests that 5-CoV RNA in these 
samples was present in insufficient amounts or was of low quality. 


Table 2. Summary of RT-PCR and sequencing results for the 14 5-CoV-positive samples. 


Sample RT-PCR Results Sequencing Results 
Deltacoronavirus Collection 
Date UDCoV HKU20 Next Generation Partial N-Gene 

Snow goose p x 
coronavirus/USA/Arkansas0009/2015 1/29/2015 N/A NA N/A 
Snow goose i 
coronavirus/USA/Arkansas0012/2015 1/30/2015 N/A N/A N/A 
Snow goose 
coronavirus/USA/Arkansas0014/2015 1/30/2015 + Ne N/A N/A 
Snow goose 
coronavirus/USA/Arkansas0017/2015 1/30/2015 B N/A N/A N/A 
Red-tailed hawk 
coronavirus/USA/Ohio1248/2015 5/3/2015 = N/A NA NA 
Mallard Did not yield CoV 
coronavirus/USA/Ohio4381/2015 8/4/2015 7 sequences N/A 
Norther shoveler 
coronavirus/USA/Mississippi8042/2015 12/22/2015 i N/A N/A N/A 

; Most closely 
Blue-winged teal 

: ae 12/22/2015 + + N/A related to 
coronavirus/USA/Mllinois2662/2017 HKU20 
Blue-winged teal 
coronavirus/USA/Missouri3057/2017 10/16/2017 H t N/A N/A 
Blue-winged teal 10/23/2017 Pi F Did not yield CoV ir ones 
coronavirus/USA/Missouri3230/2017 sequences HKU20 
Environmental 
coronavirus/USA/Maryland3464/2017 10/23/2017 ʻi 7 N/A N/A 
Blue-winged teal P 
coronavirus/USA/Tllinois2537/2017 10/16/2017 ü N/A N/A 
Yielded 85 5-CoV 
Blue-winged teal sequences from Most closely 
nae ia 10/16/2017 + + ORFla/b, S, M, related to 
coronavirus/USA/Ilinois2562/2017 N/NS7a: most closely HKU20 
related to HKU20 

Blue-winged teal 10/16/2017 $ F Did not yield CoV a a 
coronavirus/USA/Illinois2615/2017 sequences HKU20 


* N/A—Not available: was not analyzed because there was not enough RNA. 


3.4. HKU20-Specific RT-PCR and Partial N Gene Sequence and Phylogenetic Analyses of 6-CoVs 


Based on the NGS results, we designed HKU20-specific primers targeting a fragment of 
the N gene and re-screened the 5-CoV-positive samples for which we had sufficient amounts of 
RNA. Six out of 8 re-screened samples were positive in RT-PCR using HKU20-specific primers. 
Three samples yielded PCR fragments of sufficient quantity/quality and were used for capillary 
sequencing. Phylogenetic analysis of these fragments demonstrated that similar to blue-winged teal 
coronavirus/USA/IIlinois2562/2017, the three 5-CoVs were most closely related to wigeon coronavirus 
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HKU20 (Figure 3) sharing 83.38%-84.39% nucleotide identity in the N gene. The two samples positive 
in RT-PCR with UDCoV primers, but negative with HKU20 primers (Table 2) could contain insufficient 
amounts of 5-CoV RNA to generate a longer amplicon (targeted by the HKU20-specific primers) 
or could possess genetic characteristics distinct from HKU20 5-CoV and other 5-CoVs identified in 
this study. 
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Figure 3. Phylogenetic trees based on partial sequences (generated using HKU20-specific 
primers) of the N gene of blue-winged teal coronavirus/USA/Illinois2615/2017; blue-winged teal 
coronavirus/USA/Illinois2662/2017 and blue-winged teal coronavirus/USA/Missouri3230/2017. 


4. Discussion 


While y-CoVs are endemic and highly prevalent in wild and domestic birds [7-13], 5-CoVs have 
only been discovered in the US only recently. Thus, the role of wild avian species in the ecology of 
5-CoVs in the US is unknown. Environmental sampling of high-risk and co-mingling sites in Alberta 
and Saskatchewan Canada, identified 5-CoV in sparrows, that are closer phylogenetically to porcine 
5-CoVs than to those from waterfowl [19]. To date, the presence of 5-CoV is confirmed in sparrows 
in Canada and the US [13], in quail in Brazil [20] and in various birds in Australia [6], suggesting 
that 5-CoVs circulate in different avian species in the Americas. Since porcine 5-CoV often results in 
severe clinical disease and mortality in piglets, which impacts the swine industry, it is important to 
understand the morbidity and interspecies transmission rates between birds and pigs [4,5]. Due to 
their flocking behavior and abilities to fly long distances, birds can play a role in the dissemination of 
5-CoVs among themselves and other animals [12]. 

The apparent prevalence of 5-CoV in our study was only 1.13%. However, this is higher 
compared with that observed in our previous study (0.5%) that analyzed Ohio samples from the 


Viruses 2019, 11, 897 8 of 11 


period 2013-2014 [16]. Also, of these 14 positive samples, only up to 4 originated from the same 
state, suggesting that 5-CoVs spread with low efficiency in the avian species tested. Additionally, our 
inability to recover CoV sequences by NGS suggests that they are present in the samples at lower 
frequency compared with other microorganisms. The latter is more suggestive of asymptomatic carrier 
status as opposed to acute infection associated with clinical disease. These findings suggest that 
avian species likely represent a natural reservoir of -CoVs, while pigs and other mammals serve as 
spillover hosts. 

The higher prevalence of y-CoVs in the US of 4.99% is consistent with findings from some previous 
studies. In a screening study from Brazil, quail were susceptible to both ô- and y-CoV [20]. The higher 
prevalence of y-versus-5-CoVs and variations in preferred avian host species have been also reported 
by other researchers in other countries, such as Australia [6], Asia (Hong Kong and Cambodia) [11] 
and Finland [14]. This may indicate 5-CoV emerging status and its ongoing spread in the US. 

In this study, a higher prevalence of 5- and y-CoVs in aquatic vs. terrestrial birds was evident, 
with 5- and y-CoVs prevalence in aquatic birds being 1.34% and 6.3%, respectively, compared with 
only 0.6% and 0% in terrestrial birds. This suggests that aquatic birds may represent a natural reservoir 
for CoVs of terrestrial birds and pigs, and their concentration or survival may be increased in water 
sources compared with other avian habitats. Similarly, in a surveillance study in Hong Kong and 
Cambodia, 5-CoV was found in different wild aquatic bird species including gray herons, pond herons, 
great cormorants, black-faced spoonbills, and several duck species, whereas y-CoVs were found in 
little whistling ducks, tufted ducks, common teals, northern shovelers, eurasian wigeons, and northern 
pintails [11]. In wild Australian birds, ’-CoV was detected in pacific black ducks, but it was also 
detected in terrestrial birds such as curlew sandpiper, red-necked stints, ruddy turnstones, and pied 
herons [6]. This suggests that wild birds are major reservoirs of a wide range of 5- and y-CoVs, and 
the circulation of CoVs without association with clinical disease is more common than previously 
recognized [11]. However, it is important to note that 5-CoVs identified in terrestrial birds and pigs 
are more similar each other phylogenetically [12,13,21], and those from aquatic/wading birds are 
genetically more distinct. It is of interest, that prevalence of both y- (92%) and 5-CoVs (86%, Table 1) 
was higher in colder (October—February) than in warmer (March-September) months. This indicates 
that, similar to previous findings on avian, animal and human CoVs, there are may be seasonal 
(and migration associated) fluctuations in the prevalence of avian coronaviruses in wild bird in the 
US [22-25]. 

Because phylogenetic analysis of the partial N gene of the 4 newly identified 5-CoVs (Table 2) 
revealed that they were most closely related to HKU20 (wigeon) 5-CoV, we hypothesized that it may 
be the parental strain recently introduced into the US. 

In our study, these newly identified 5-CoVs were more closely related to other 5-CoV from aquatic 
birds, but not terrestrial avian species and pigs, which can be explained by the fact that aquatic birds 
occupy a separate ecological niche, while terrestrial birds and pigs may share some co-mingle sites. 
Although, we did not identify 5-CoV strains genetically similar to sparrow and porcine 5-CoVs, it is 
important to note that porcine 5-CoV outbreaks in the US were predominantly detected in the states 
with the highest density of pig population that in turn have a significant geographical overlap with 
the Mississippi Flyway (Figure 4) [26]. While a small pilot study failed to identify porcine 5-CoVs in 
wild waterfowl in the Mississippi Flyway [27], the overlap of the bird migratory pathways and high 
density swine farms creates favorable conditions for birds and pigs to exchange 5-CoV pools directly 
or through some other intermediate hosts. 
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Figure 4. Porcine deltacoronavirus was detected in at least 18 US states (shaded in blue) along or 
adjacent to the Mississippi Flyway (blue contour). Orange circles mark states with the highest pig 
inventory [28]. Black duck icons identify the states where avian 5-CoVs were found (this study). 


5. Conclusions 


In conclusion, this is the first report of the presence of HKU20-like 5-CoVs in different avian 
species in the US. The close relatedness of all the strains we were able to sequence to HKU20 suggests a 
recent introduction of 5-CoVs in the US and identification of this strain as potentially parental. We also 
demonstrated that aquatic birds are infected with 5-CoVs and y-CoVs more frequently than terrestrial 
avian species. This, together with the observation that porcine 5-CoVs are more closely related to 
5-CoVs identified in terrestrial birds, suggests that waterfowl might represent a natural reservoir 
for 5-CoVs. Although y-CoVs were more prevalent than 5-CoVs, consistent with previous studies, 
both CoVs were detected more frequently during the cold season in our study. Further studies are 
necessary to investigate the role of aquatic bird 5-CoVs in the epidemiology of 5-CoVs in swine and 
terrestrial birds. It remains to be established: (i) if/how long avian 5-CoVs were present in South and 
North America prior to the porcine 5-CoV outbreaks in the period 2013-2014; (ii) which avian species 
represent significant natural 5-CoV reservoirs; and (iii) what is the potential of avian 5-CoVs to cross 
the interspecies barrier and infect swine. 
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